

% ratio of Px/Py vs degree of anisotropy simulation data
% 1:1, 2:1, 3:1, ... , 10:1 degree of anisotropy
% 2.2 um diameter - 10nm radi force face - 45nm thickness ... 
% 200nN force - Px1000 & Py1000 
y=[1 1.229 1.392 1.525 1.635 1.728 1.810 1.879 1.934 1.982]; %Px/Py values for each DoA  
x=[0 1 2 3 4 5 6 7 8 9]; %DoA
polyan = polyfit(x,y,4); %Fitting a polynomial for data
xan=0:0.001:9;
yan=polyval(polyan,xan);

% plot(x,y);
% hold on;
plot(xan,yan);
hold on


plot(7.802,1.927, 'r*'); %plotting red asterix




